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The radiative muon capture (RMC) on a proton is analyzed by means of 
heavy baryon chiral perturbation theory. The emitted photon energy spectrum 
is calculated and compared with the experimental data by taking the spin 
sum on the muonic atom states. We find that one-loop order corrections 
■ to the tree order amplitude modify the photon spectrum by less than five 

percent. This calculation supports that the theory is under a quantitative 
control as far as the chiral perturbation expansion is concerned and indicates 
that the discrepancy between the pseudo-scalar coupling constant required 
by the RMC experiment and the one deduced from ordinary muon capture, 
^ . the value of which is also supported by chiral perturbation calculations, will 

>- ! remain unexplained from the theoretical side. 
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The induced pseudo-scalar coupling constant gp was determined from the ordinary 
muon capture (OMC) reaction on a proton (p + // — > n + i/)|TJ to be gp = 8.7 ± 1.9. 
Despite its large error bar, this value is clearly consistent with the theoretical prediction 
by Bernard et al. M using heavy baryon chiral perturbation theory (HBChPT), gp(q 2 = 
—0.88m 2 ) = 8.44 ± 0.16. This is also comparable to the value that Fearing et al. @ 
evaluated by means of HBChPT in their work on OMC, # P (-0.88m;;) = 8.21 ± 0.09. 
Together with the PCAC prediction gp CAC {— 0.88m 2 ) = 8.42, all theoretical investigations 
agree with the experimental result on OMC. However, the momentum transfer involved 
in OMC is far from the pion pole, q 2 = m 2 , where the pseudo-scalar coupling should play 
important role in the reaction amplitude, accounting for the large error bar. 

Since RMC involves a momentum transfer closer to the pion mass, it is considered to 
be more suitable to measure the constant gp and lower the error bar. For this purpose, the 
photon energy spectrum from the radiative muon capture (RMC) (p + fi — > n + u + j) has 
been measured in TRIUMF [[| and compared to the model prediction of ||. Surprisingly, 
the experimentally detected photon spectrum could be explained only if the pseudo-scalar 
coupling constant is enhanced in the model by a factor of 1.5 relative to the value given 
by PCAC H or that determined in OMC. A calculation to tree order, recently reported 
by Meissner et a/.0, further confirms this discrepancy. 

The purpose of this paper is to see whether or not this discrepancy can be eliminated 
by higher order terms in the treatment of the strong interaction sector of the process. 
It is natural to ask whether any important Feynman diagrams have been ignored in the 
phenomenological model, in particular in light of the direct chiral perturbation calculation 
of the pseudo-scalar coupling constant by Bernard et al. [[|] which agrees with the PCAC 
prediction. Experiments are currently being planned || to increase the precision. 

In this work, we shall calculate the RMC amplitude and the photon energy spectrum 
using the HBChPT up to the next-to-next to the leading order (N 2 LO), that is, to one- loop 
order. We shall also investigate the photon energy spectrum by taking various ansatze on 
the spin states of the muonic atom. 

Heavy baryon chiral perturbation theory(HBChPT) H] provides a systematic way of 
making, in the presence of nucleons, a chiral perturbation expansion in powers of Q/A x , 
where Q is a typical momentum scale and/or the pion mass and A x is the chiral symmetry 
scale, A x ~ lGeV. Since the momentum transfer involved in RMC can be of the order 
of the muon mass, m M ~ 0.1A X , the chiral expansion is expected to converge sufficiently 
rapidly. 

The Feynman graphs contributing to RMC can be classified into two classes as shown 
in Fig. 1: (a) the first corresponds to those graphs where the photon line is attached to 
the lepton, therefore, leaving the nucleon line to form the 3-point vertex of WNN (weak 
current-nucleon-nucleon) , (b) the second corresponds to the graphs where the photon is 
attached to the nucleon line and the vertex with the exchanged pion coupled to the weak 
boson, which is schematically a 4-point vertex of VWNN (electro-magnetic current-weak 
current-nucleon-nucleon). Indeed, those two vertex graphs shown as blobs in Fig. 1 can 
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(b) 
(VWNN) 



Figure 1: Diagrams for three and four point Greens functions. 



be expressed in terms of Green's functions as follows, 

/ '''''' 



3 r d 4 x- 

Jp(Qj> k ) = -XnSHrHtf) II J ?2^ exp ' Xl ~ k ' ' X2 + k ' X3 ^ 



(2nY 

x(0|T W^x 1 )N(x 2 )N(x 3 )\0)S N 1 (k) Xp , (1) 
4 r d 4 x- 

M* f p(q, qw, k) = -ix{Sj; l {k') \{ J ^^ ex P H(~9 ■ x 1 + q w ■ x 2 - k' ■ x 3 + k ■ x A )] 
x(0|T V:(x 1 )W^x 2 )N(x 3 )N(x,)\0)S N \k) Xp , (2) 



where e,f and a, (3 are the iso-spin and Lorentz indices, respectively, and T stands for the 
time ordering on the currents and fields appearing on its right. Note that the vacuum 
expectation value in the above equations is the Green's function. The four-momenta 
carried by proton, neutron, neutrino, muon are denoted by the particle symbols, p, n, 
v, fi, respectively, whereas the four-momentum of the photon by q. The momentum 
transfers described in Fig. 1 can be written as qw = H — v and qj = qw ~ 1- Here Xp 
(Xn) is a two-component spinor of the proton (neutron) with the normalization condition^] 
X^X — E + rriN, and S^ f 1 (k) is the inverse of the heavy nucleon propagator to be derived 
below. In the above, k (k') denotes the residual momentum of the proton (neutron) so 
that 

p» = m N v^ + k^, (3) 
n" = m N v^ + k'^, (4) 

with the velocity four- vector = (1,0). is the physical nucleon mass. The heavy 
nucleon field N is defined from the nucleon field ^ n as 

N = e imNV - x P + m N) (5) 

where P + is the projection operator defined by P + = ~(1 + v • 7). 
3 Actually x depends on ?A Here we have set i/' = (1,0). 
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The Green's functions in Eqs.(|T]) and ([|) are to be obtained from the effective chiral 
lagrangian with nucleons and pions, the expression of which reads 



C = C + A + C 2 + 



(6) 



where Co is the leading order lagrangian given in JTof and C\ is of the 1/rajv correction 
(NLO) which will be specified below. C 2 is the next-to-next leading order (N 2 LO) effective 
lagrangian. The ellipses stand for higher order lagrangians irrelevant for our calculation. 
C\ reads [|i~0| 



1 



-N 



2m N 
+b 1 Tr(x+) + b 2 (x 



-D 2 + (v ■ Df + 2g A {v ■ A, S ■ D} 



S v ] (2(1 + h)A„A„ + z(l + 6 6 )/+ + i(l + b 7 )v 



N, 



(7) 



where we have adopted the notations D M , A M , x+, and i>£, as defined in Ref.flO 
except for an additional multiplication factor 1/2 for our /+. g A is the axial- vector 
current coupling constant, and bi are the low energy constants that cannot be fixed by 
the theory, but will be determined from experiments. We should note that for RMC, 
only the two constants b^ and bj are relevant. They are determined from the anomalous 
magnetic moments of nucleon as b§ = = 3.71 [], 67 = ks = —0.12, where Ky and ks 
are the iso-vector and iso-scalar anomalous magnetic moment, respectively. The N 2 LO 
lagrangian C2 containing low energy constants and an anomaly term reads 



lot 



(2) 



T^VVWC/t + R^V U U^V U U) + a$Tr (L^UR^lP) + C wz 



4:171% 



N 



iD a v ■ DD a - 2i[S a , S p ]D a v ■ DD 8 + ig A [2D ■ AS ■ D - 2D a S ■ AD Q 



+2S ■ DA ■ D + 2v ■ DS ■ Av ■ D - 2{S ■ D, {v ■ D, v ■ A}} 



abed 



D a A b D c v d ] 



+i[-v a D p + v p D a - 2v a [S 13 , S-D} + 2v?[S a , S ■ D] + 2g A (v a S p - v p S a )v ■ A)B af3 
+iB aP [v a D p - v p D a - 2v a [S p , S-D} + 2v f} [S a , S ■ D] - 2g A (v a S p - v p S a )v ■ A]] N 



c 3 v a [D^, f+] + c 4 [S a , S^]{v ■ D, f+ } + c 8 g A v a ' 



(47TA) 

+c 12 [S a , S p ]{v ■ D, vlg} + c 13 g A S a [D a , f~J + ic lA g A S a [D a , x 



N, 



with 



.h ,+ .67 s 

1 2 J a/3 + 1 2 V a/3 



(9) 



where Cwz is the Wess-Zumino lagrangian [|i~4"| . Note that we have eight low-energy con- 



stants, among which seven of them are needed for this work: a^ 1 + a\o = 1.43 x 10 



.,(2) 



1-3 



'Up to N 2 LO, Ky is renormalizcd by b$ — 



Ky- 
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determined from a rare pion process [11], C3 = 5.34 (C13 = 2.37) from the iso-vector vec- 
tor (axial-vector) radius, cu = —1.37 from the Goldberger-Treiman discrepancy^, |3|, 



c 8 
P 



3.27 from the A contribution to En? ||12||, C4 = —22.27 and cyi = —0.79 from the 



uj and A contributions to P,' [13 



In short, the constants are completely determined 
for calculating up to the N 2 LO order in chiral perturbation expansion. 

We are now in a position to calculate the relevant Feynman graphs. The chiral power 
counting rule for A-nucleon processes is that for a Feynman graph with Vi vertices of 
type i, L loops, and C separately connected pieces, the power index of Q is given by 
rj = 4 — A — 2C + 2L + YU Vi^i with Aj = di + rij/2 — 2, where rii is the number of 
nucleon lines and di is the number of derivatives or powers of m v at the i— type vertex. 
In the presence of an external gauge field, Aj is constrained by chiral symmetry to be 
Aj > — 1|TIJ. Thus the leading order of matrix elements J is 0(1) and that of M is 
0{Q^ 1 ). However, the leading order amplitudes of Fig. 1(a) and Fig. 1(b) are of the 
same chiral order, because the muon propagator in Fig. 1(a) is of order Q _1 since it 
carries the photon momentum in the denominator. 

Now we split the weak-current into the V — A form so that the J and M can be written 



J 

M 



Jv — Ja, 
M v - M A , 



The most general forms of J turn out to be, (with = (1, 0)) 

J v (qj, k = 0) = ft, Jyfa, k = 0) = iax qjfX + qjf, 
J\(qj, k = 0) = a ■ qjff, J A {qj, k = 0) = aff + qjd ■ qjf. 



A 
2 ' 



(10) 

(11) 



(12) 



where as is done in what follows, the initial and final state nucleon spinors are omitted. 
ff and ff denote the nucleon vector and axial-vector form factors, respectively. They 
read as, up to N 2 LO, 



f? 



1 + 



c 3 



(4vrA 



18(47T/ ff ) 2 ^ (47Tf n 



2 ~(l + 2gl)ml- 1 ^ 
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Am 2 



N 



- + K V )qj, 



fl 



2m 7v 
1 

2l7lN 



[1 + K V )\qj\ + 



9a 



Qj\Qj\ + 



9a \ n(4m 2 n -qj) 



4m„ 



™o(Qj)\Qj\ 



\QJ\ 



9a 

9a 
9a 



1 + 



Cl3 



Q 2 j 



N , 



Cl3 



A, 



2m 



[1 + K{qj)q 2 j]\qj\ 



1 - 



2m 2 r ci4 




\Qj\ 



N 



fo(qj) 

(13) 

(14) 
(15) 
(16) 

(17) 
(18) 
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with 



1 2 

/„(?) = / dx\u[l-x{l-x)^l (19) 

ml 



f 1 1 

7710(g) = 1 - / dx (20) 
"'o Wl + x(l — x)q 2 /ml. 



The common factor 2m^ is omitted in the above equations, and qa = 1.25. The conver- 
gence of the form factors is found to be quite good as is discussed in Ref. @, ||, [10]. 



Under the Coulomb gauge, the renormalized inverse propagators for our calculation 
can be written simply as 

S N \k) = v .k + -^—[k 2 -vk% (21) 
A^(g) = q 2 -m\. (22) 

We choose the coordinate frame such that the neutrino lies in the z-direction and the 
photon in the x-z plane, respectively, i.e., v = (0,0, 1) and q = (sin#, 0, cos9), where 6 is 
the angle between the neutrino and the photon. Then we can decompose M into so-called 
reduced amplitudes for each muon spin states lying along z-axis, m s = ±1/2 JT6|j , 

M = -e T ■ M T 5 msil/2 + e Li3 M L 5 ms ,_i/ 2 , (23) 

where = ef 1/2) and e{ = ef_ 1/2) with ef ma) = u u ^ 5 {l - ^)u^ s \ 

Then generally one can decompose them into different spin operators as specified in 
Tables 1 and 2, 

6 10 

-E T ■ Ma,T = E A i°a,i, M A,L = E A i°a,i, ( 2 4) 
i=l i=7 
9 13 

-e T ■ M V)T = E M y , L = E BiO b4 , (25) 



i=l i=10 



where O a ^ and Ai (Ob,i and Bi) are operators and corresponding form factors, respectively. 
The Ward-Takahashi (WT) identity was found to be quite useful in reducing redundancy 
in the form factors. Summation runs over all possible effective operators. 

Some LO and NLO contributions in Ai in Table 1 contain the pion propagator taken 
at qw, which make the difference between RMC and OMC We present the results in Tab. 
1 and 2 calculated in Coulomb gauge. The formulae for the matrix elements are quite 
lengthy and uninstructive; we leave their explicit expressions to a forthcoming paper fllT 



For the contribution of N 2 LO we give the values of their maximum among the entire 
range of photon energies. Among the contributions of this order, the most important 
one comes from the intermediate excitation of a A contributing to the term proportional 
to C4. We have multiplied by m M in the last column of Table 1 to make the numbers 
dimensionless. One can see that at their maximum, some of them are comparable with 



6 




Figure 2: Photon spectrum of contributions of each order. 
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Figure 3: Photon spectrum for different spin states of muonic atom. 
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O b ,i 


Af° 


a ATT n 


i /i Loot) \ 
l^fi-^-i \max 


1 


iq ■ e* x e T 







0.030 


2 


a ■ e* €t ■ q 





-2^(! + ^) 


0.207 


3 


a ■ q e* ■ Or 


-g A A^{q J )E 1 


3^(1 + 


0.232 


4 


a ■ q w e* ■ e T 


g A A„{qj)E v 





0.001 


5 


a-q e* -q w e T ■ q 








0.001 


6 


a-q w t -q w e T -q 








0.001 


7 


ie* • q x q w 





x(l + A^{q w )m^E v ) 


0.033 


8 


— * 

a ■ e 


-5f j4 A 7r (g l y)m^ 


[(! + "*) 
x(l + A n (q w )m^E u )y 

-(1 - A 7r (g H /)m M E 7 )] 


0.219 


9 


a-q? - qw 


g A A n (qj)E^(l+ 
2A 1T (q w )m tl E v ) 


x(l + A^{qj)m^E v ) 


0.184 


10 


a ■ q w e* • q w 


-g A A 7T (q J )E I/ (l+ 
2A 1T (q w )m tl E v ) 





0.015 



Table 1: Operators and form factor Aj for the M A . where E 1 (E u ) is the photon (neutrino) 
energy. And y = cos#, 9 is the angle between the neutrino and photon. 



the 1/rriN corrections. However, in the total spectrum, their correction is less than five 
percent as is shown in Fig. 2. 

We are now ready to discuss the results of our work. We start by looking at the 
role of the pion propagators. The momentum transfer, qj, is always space-like, i.e., 
q 2 j ~ — qj, due to the on-shellness of the incoming and outgoing nucleons. This is the 
reason why A n (qj) is suppressed with the important contribution coming from f{ and 
fi instead of from f£ and f£ in Eqs. C[l3| , [T6| , [T7| , p~8D . On the other hand, q^ increases 
almost linearly with E 1 and becomes time-like when E 1 is greater than ~ 50 MeV, since 
q^r ~ 2m ii E 1 — m 2 . Note that this is the region where the photon energy spectrum is 
established in the experiment. Hence A v (q\y) is enhanced in a high photon energy region 
while A n (qj) is always suppressed. Consequently the LO distribution comes mainly from 
the three terms fi' LC \ ff' LO and A%° . In particular, the contribution from Aq°, the 
so-called Kroll-Ruderman (KR) term, carries about thirty five to sixty percents of the 
photon spectrum for E 1 > 60MeV. 

In Fig. 2 the spin averaged photon energy spectra for the LO, NLO and N 2 LO con- 
tributions are plotted. We find that the result of the phenomenological model || can 
be more or less reproduced by the LO and NLO contributions. For the experimentally 
measured region of the photon energy, the NLO contribution remains within 20 % of the 
LO contribution. 

To summarize, we found that the next-order correction to the NLO description is 
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i 




j^NLO 


\'"'H ±J i \max 


1 




1 

2m /v 


0.080 


2 


t -q w e T -q 





~ 10" 4 


3 


ia ■ (e* x e T ) 




0.007 


4 


ia ■ (qx q w ) e* • e T 





0.003 


5 


• (e* x gV) e T • q 





0.003 


6 


ia ■ (e* X q) or ■ q 





0.009 


7 


-ia ■ (e T x g) e* • g w 





0.003 


8 


-ia • (e T x q N ) e* • g^ 





0.010 


9 


ia - (qx q w )e* ■ q w e T ■ q 





~ 10" 4 


10 


e* ■ g w 


1 

2m /v 


0.029 


11 


ia ■ (qx q w ) e* • g 





0.020 


12 


• e* x giy 


2m„( 1 + ^) 


0.014 


13 


i? • 6* x g 




0.010 



Table 2: Operators and form factor A\ for the My: Operators are identical to those 
occurred in Ref. |p~6[| 



negligible and does not remove the discrepancy present at that order: The further correc- 
tion to the spectrum does not change appreciably the results of the previous theoretical 
calculations |)]. [|. 

This may seem disappointing in the sense that the persistent puzzle is not resolved by 
our higher order calculation. On the other hand, our calculation is tightly under control 
and the fact that the next-order terms to the NLO contribution are negligible implies that 
our theoretical treatment has converged. It is then legitimate to ask what mechanisms 
other than strong-interaction dynamics could be the cause of the discrepancy. As an il- 
lustration of such alternative mechanisms, we have considered the effects of various spin 
states in which the muonic atom could be formed. The possible photon spectra for these 
spin states are given in Figure 3. It is interesting to see that if one assumed only the triplet 
state of the atom to be occupied, then one would reproduce the observed photon energy 
spectrum. While we are not claiming that this could account for the discrepancy, such 
non-strong interaction mechanisms could not be ruled out. Given the theoretical confi- 
dence in calculating higher-order chiral corrections, it seems imperative that the presently 
available experiment be re-scrutinized or that more refined measurements be made before 
concluding that the constant gp is so drastically deviating from the Goldberger-Treiman 
value. 

This work is supported in part by the Korea Science and Engineering Foundation 
through Center for Theoretical Physics of Seoul National University, and in part by Korea 
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